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ABSTRACT 

The beginning of photoionization marks the transition between the post-Asymptotic Giant Branch 
(post-AGB) and planetary nebula (PN) phases of stars with masses < 8 Mq. This critical phase is 
difficult to observe, as it lasts only a few decades. The combination of jets and magnetic fields, the 
key agents of PNe shaping, could give rise to synchrotron emission, but this has never been observed 
before in any PNe, since free-free emission from the ionized gas is expected to dominate its radio 
spectrum. In this paper we report radio continuum observations taken with the Australia Telescope 
Compact Array between I and 46 GHz of the young PN IRAS 15103—5754. Our observations in 2010- 
2011 show non-thermal emission compatible with synchrotron emission from electrons accelerated at 
a shock with spectral index a ~ —0.54. However, in 2012, the spectral index a ~ —0.28 is no longer 
compatible with synchrotron emission in these types of processes. Several hypothesis are discussed to 
explain this change. The more plausible ones are related to the presence of the newly photoionized 
region in this young PN: either energy loss of electrons due to Coulomb collisions with the plasma, or 
selective suppression of synchrotron radiation due to the Razin effect. We postulate that the observed 
flattening of non-thermal radio spectra could be a hallmark identifying the beginning of the PN phase. 


1. INTRODUCTION 

Planetary Nebulae (PNe) represent one of the final 
phases in the evolution of stars with masses < 8 Mq, 
when they become hot enough to photoionize the cir- 
cumstellar envelope ejected during the Asymptotic Giant 
Branch (AGB). This nebular photoio nization proceeds 
very quickly, over only a few decades (|Bobrowskv et ahl 
[19^ . making it extremely difficult to study the events 
produced just at the transition to the PN phase, which 
are crucial to determine how PNe form. 

Many PNe exhibit complex axisymmetric or multi¬ 
polar morpholo gies that are believed to be caused by 
collimated jets (|Sahai fc Trauged 1 199811 launched during 
the post-AGB (which just precedes the photoionization 
of the envelope) and PN phases. The exact mecha¬ 
nism for the launching of these jets is still unknown, 
but p ossible scenarios are the presence of magnetic 
fields (IGarcfa-Segura et aO2005ll and/or binary systems 

(iSoker fc RaDDaDortll2000riNordhaus &: Blackman! 120061 : 

Ide Marcoll2009ll . The combination of jets and magnetic 
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fields lead to the presence of synchrotron radiation in 
seve ral astrophysical conte xts, such as acti ve Kalo-ctic nu- 
clef (Bridle fc PerI^ll984D . mic roouasars (iMirabel et ^ 


), or massive protostars (|Carrasco-Gonzalez et al 


2010f) . This emission is caused by the motion within a 


magnetic field of electrons accelerate d to relativistic ve - 
locities at shock fronts induced by jets (IAchterbersj2000[l . 

Magnetic fields have been det ected in the c i rcum- 
stellar envelopes of some PNe (iMiranda et al.1 120011 : 
iSabin et al.ll2007ll . but synchrotron emission has never 
been unambiguously reported in PNe before (see, e.g., 
ICasassus et al.ll200A for a review of previous attempts 
for such a detection), even though the energy condi¬ 
tions in PNe are w ell-suited for the production of syn- 
chroH on radiation (iDgani fc Sokei][r998t ICasassus et al.1 
l2007f) . The reason is that the strong free-free radiation 
from ionized gas that characterize these objects is ex- 
pected to overwhelm any type of non-thermal emission. 
ICasassus et al.1 (1200711 claimed that synchrotron emission 
is a possible source of centimeter-wave excess in the ra¬ 
dio continuum flux of some PNe, over what would be 
expected from free-free emission alone, but the avail¬ 
able data could not confirm this hypothesis. On the 
other hand, non-thermal emission has been observed in 
a few post-AGB sources fcohen et'aril2006l: iBains et'aLl 
120091 : IPerez-Sanchez et al.ll2013ll . where photoionization 
has not yet started but jets are already present. All this 
indicates that the temporal window for a possible study 
of synchrotron radiation in PNe is extremely short, since 
it cannot be detected very shortly after the onset of pho¬ 
toionization. 

The source IRAS 15103—5457 (hereafter 115103) is an 


optically obscured PN showing a high-velocity jet traced 
by w ater masers (|Bendiova et al. 1 120141 : iGomez et al.1 
l2015f) . Water maser jets are seen in a few late AGB 
and post-AGB st ars, and are termed “water fountains” 
(ITmai et aLl[2(iol . 115103 is the first known case of a 
“water fountain” that has already entered the PN phase. 
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Table 1 

Observations performed 


Date 

Frequency®' 

(GHz) 

Array'^ 

2010-0ct-07 

5.5, 9.0 

H214 


43.0, 45.0 

H214 

2011-Aug-ll 

21.8, 24.0 

H168 

2012-Oct-23 

2.1 

H214 


5.5, 9.0 

H214 


22.0, 24.0 

H214 

2012-Dec-10 

2.1 

6B 


5.5, 9.0 

6B 


26.0 

6B 


^ Central frequency of each band. 
^ ATCA configuration used. 


suggesting that it may be one of the youngest PN known. 
Therefore, it is an excellent candidate to study the pro¬ 
cesses taking place during the transition into this phase. 
Moreover, interferometric radio continu um observations 
obtai ned in the Red MSX Survey (RMS, iLumsden et al.l 
l201.‘lfl show a flux den sity lower at 3.6 cm than at 6 cm 
(jUraiihart et al.l[2007[l . which suggests the emission may 
be non-thermal. A confirmation of this trend was needed, 
since only two wavelengths were observed in this sur¬ 
vey, and problems such as calibration uncertainties or 
extended flux missed by the interferometer may have af¬ 
fected these values of flux density. 

In order to confirm the presence of non thermal emis¬ 
sion in 115103, we present new radio continuum obser¬ 
vations covering a much wider range of wavelengths in 
different epochs. In Section 2 we describe the obser¬ 
vations performed. In Section 3 we explain the results 
derived from the observations. Section 4 is devoted to 
the different hypothesis we analyse to explain the varia¬ 
tions observed in the spectral energy distribution (SED). 
These hypotheses are put in context and the predictions 
of the future evolution of the SED are presented. Section 
5, finally, presents the conclusions of this work. 

2. OBSERVATIONS AND DATA REDUCTION 

Radio observations towards 115103 were performed us¬ 
ing the Australia Telescope Compact Array (ATCA) in 
several epochs between 2010 and 2012. Table [T] sum¬ 
marizes the different frequencies and array configuration 
used. All observations used the Compact Array Broad¬ 
band Backend (CABB), which provides two independent 
intermediate frequency (IE) outputs of 2 GHz in dual lin¬ 
ear polarization. In the particular case of observations 
around 2.1 GHz, both IE were centered at the same fre¬ 
quency, since each IE covers the whole frequency band of 
the receiver (to provide redundancy in case of correlator 
failure), so we only processed one of them. The obser¬ 
vations on August 2011 used the 64M mode of CABB, 
which samples each 2 GHz bandpass into 32 frequency 
channels of 64 MHz. All other observations used the IM 
mode, which samples 2048 channels, each of them 1 MHz 
wide. 

The complex gain calibrator was PKS 1613-586 in all 
datasets. The absolute flux calibrator was PKS 1934-638 
at all observed frequencies, except at the beginning of the 
observations centered at 5.5 and 9.0 GHz on December 


2012, when this source was not yet above the horizon, 
and PKS 0823-500 was used instead for this purpose. 
Since PKS 1934-638 is the recommended source to set 
the flux density scale in ATCA data, in these particular 
set of data, we independently calibrated the data that 
used PKS 1934-638 and PKS 0823-500 as flux calibrators, 
and rescaled the latter data (by a factor of ~ 1.10— 1.15, 
depending on frequency) to achieve the same final flux 
density on the complex gain calibrator throughout the 
observations that day. This indicates that the flux den¬ 
sity model of PKS 0823-500 used by the data reduction 
package was in error by ~ 10 — 15%, as pointed out in 
the ATCA calibrator catalog. The two sets of uv-data 
were merged after this rescaling. The flux calibrators 
were also used to calibrate the bandpass at frequencies 
below 15 GHz. At higher frequencies, we initially used 
PKS 1921-293 to correct the bandpass, but checked the 
results against PKS 1934-638, to confirm that the result¬ 
ing spectral slopes in our data were correct. 

All data were calibrated and processed with the 
MIRIAD package using standard procedures. Calibra¬ 
tion, initial imaging and self-calibration were performed 
across the entire 2 GHz bands. Final uv-data were then 
imaged in successive 200-MHz intervals to obtain the 
intra-band variation of flux density with frequency. The 
two exceptions for this processing mode were the data 
between 43—45 GHz, where the signal-to-noise ratio was 
lower and we present the average of each 2 GHz band, 
and the data between 21-25 GHz taken on August 2011 
where we kept the original 64 MHz spacing provided by 
the correlator. To double check our results, we repro¬ 
cessed the data by splitting the raw uv-data into the 
same frequency intervals, which were then independently 
calibrated, self-calibrated, and imaged. The resulting de¬ 
pendence of flux density with frequency was the same 
using both methods of splitting. In all cases, broadband 
data were weighted with robust parameter 0.5 and fre¬ 
quency synthesis was applied, using MIRIAD task “in¬ 
vert” . Images were subsequently deconvolved with the 
CLEAN algorithm, as implemented in task “mfclean”. 

Low-frequency observations are contaminated by ex¬ 
tended Galactic background emission. This affects the 
shortest baselines. In order to obtain reliable maps we 
had to restrict the data to baselines > 6 — 8 kA at fre¬ 
quencies < 3.5 GHz. With this restriction, the 2.1 GHz 
data taken in October 2012 with the compact H214 ar¬ 
ray were unusable. We also note that a visual inspection 
of the maps between 4.5 and 6.5 GHz obtained on Octo¬ 
ber 2010, revealed significant artifacts and a higher noise 
(a factor of 3) than those at > 8 GHz. We think that 
they could be affected by low-level radio-frequency inter¬ 
ferences (RFl), which we could not identify in uv-data 
for flagging, but that affects the data reliability. Fur¬ 
thermore, the corresponding flux density points show a 
relatively large scatter. 

To reconstruct the SED behavior with time we have 
used the archival and published data summarized in Ta¬ 
ble [21 These dat a have been taken with ATGA by other 
authors in 1991 (Ivan de Steene fc Pottaschlll993ll. 1999 
(iGaensler et al.l l200r)tl . 2002 and 2004 ( Urouhart et al.l 
l2007f) , as well as with the Molonglo Observatory Synthe¬ 
sis Telescope in 1998 for the Second Mol onglo Galactic 
Plane Survey (1998) ([Murphy et al.ll200^ . We used the 
published flux densities for the data from 1991 and 1998. 
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Table 2 

Archival data 


Date 

Frequency®' 

(GHz) 

Flux density 

(Jy) 

rms^ 

(Jy) 

References 

1991.4 

4.8 

0.1145 

0.0008 

1 

1998.4 

0.8 

0.193 

0.006 

2 

1999.2 

1.4 

0.199 

0.003 

3 


2.5 

0.219 

0.004 

3 

2002.7 

4.8 

0.1509 

0.0006 

4 


8.4 

0.1088 

0.0004 

4 

2004.9 

4.8 

0.0989 

0.0007 

4 


8.4 

0.068 

0.002 

4 

2012.0 

1.8 

0.052'^ 

0.002'= 

5 

References. — 1. 

Ivan de Steene ^ Pottaschl (Il993l). 2. 

IMurobv et 

al.l fl2007l). 3. IGaensler et al.l (120001) 4. lUrauhart et al.l 


^ Central frequency of the observation 
^ One-sigma rms noise of the map 

^ Mean value for the frequency intervals in which we split the whole 
band 


We note that data from 1991 do not have enough uv cov¬ 
erage to obtain images, so their flux density values may 
have a significantly larger error than quoted. We down¬ 
loaded and processed the ATCA archival data on 1999 
and 2002 to obtain the values in Table [H The fluxes 
on 2004 were measured on processed images provided 
by the authors (jUroiihart et all 120071 1 . The data taken 
in Janu ary 2012 centered at 1.8 GHz for the MAGMO 
project (|Green et al.ll201^ . included 115103 within their 
primary beam, with their phase center located 13.7' away 
from our source. Therefore, we processed these data cor¬ 
recting by the response of the ATCA primary beam as 
implemented in the MIRIAD task “linmos”. 

3. RESULTS 

In Figure [T] we present the radio continuum emission in 
115103 at different frequencies and epochs, obtained with 
our observations and the archival data. Our radio contin¬ 
uum observations show a dependency of flux density S^, 
as a function of frequency that is atypical for a PN. Ra¬ 
dio continuum emission from PNe normally consists in 
free-free (bremsstrahlung) radiation in the ionized gas, 
and follows a relation with frequency 3,^ cx z/“, where 
the spectral index a ranges from a ~ 2 in the optically 
thick regime at low frequencies (usually u < 10 GHz), to 
—0.1 for optically thin emission at high frequencies (usu- 
^ 10 GHz). The turnover frequency separating 
these two opac ity regimes is tho ught to depend on the 
age of the PN (iKwok et al.lll98lll . 

Howeyer, our data taken in the first two sets of obser- 
yations (2010-2011; black and green solid lines and black 
filled circles) show radio continuum emission with a steep 
decline at high frequencies (Figure [T|). We found a spec¬ 
tral index a = —0.54 ± 0.08 between 8 and 25 GHz, in¬ 
dicating a non-thermal origin for the emission. The data 
at higher frequencies (42-46 GHz) show an apparently 
steeper spectral index, but giyen that it is obtained from 
only two points, its significance is low (a ~ —3 ± 3). 
At lower frequencies (4-5.5 GHz), the data haye lower 
quality, with artifacts in the maps that may indicate the 
presence of low-leyel RFI, as explained aboye, so the pos¬ 
sible turnoyer suggested by the data at these frequencies 


is not reliable. 

Our most recent ATCA data (October-December 
2012), taken about one year later, show a significantly 
flatter spectrum oyer the whole range from 1 to 24 GHz. 
In particular, the data taken on October 2012 (blue solid 
lines in Fig. [T]) follow a spectral index of a = —0.28±0.08 
between 4.5 and 25 GHz. The slope is similar for data 
taken on December 2012 between 1 and 10 GHz (orange 
solid lines in Fig. [T]). 

The data taken on December 2012 allow us to resolye 
some structure in the radio continuum emission. Figure 
[2] (top) shows a map of the emission integrated between 
8 and 10 GHz. The emission is dominated by an unre- 
solyed intense core, but there is much weaker emission 
extending along three arms toward the northeast, south¬ 
west, and southeast. Note that the contour leyels In this 
map are not linear, but proportional to powers of 2, and 
the highest leyel nearly corresponds to the half power 
leyel. Figure [2] (bottom) shows a spectral index map ob¬ 
tained with the emission between 8 and 10 GHz (map in 
Fig. [2]-top), and the radio continuum emission between 
4.5 and 6.5 GHz, both conyolyed to a common angular 
resolution of 2.25". As expected in yiew of the SED, 
most of the emission is non-thermal (with a < —0.1), 
including the central bright core. Howeyer, there are 
some areas that are compatible with thermal emission 
(a > —0.1), specially toward the northeast (the posltiye 
spectral Indices toward the southeast are at lower inten¬ 
sity leyels and are less reliable). The intensity of this 
possible thermal emission is more than 50 times weaker 
than the non-thermal peak, and it does not significantly 
affect the deriyed spectral Index from Fig. [TJ 

The time eyolution of flux density deriyed from archiyal 
obseryations is also consistent with the emission becom¬ 
ing weaker with time, a yariation that seems to be dra¬ 
matic at low frequencies. The only exception to this his¬ 
torical weakening Is the data from the MAGMO project, 
taken In 2012, and optimized for spectral line work. The 
MAGMO data suggest a turnoyer, and a slightly lower 
flux density than our obseryations at these frequencies, 
which were taken later. Howeyer, we consider these 
flux densities less reliable, giyen the relatlyely large dis¬ 
tance of the source from the phase center of the obserya¬ 
tions, and the use of a different primary flux calibrator: 
PKS 0823-500 Instead of PKS 1934-638 In our data. As 
noted aboye, the flux density model of PKS 0823-500 In 
MIRIAD might be incorrect, and so the combination of 
positional offset and calibrator flux density error margin 
reduce the significance of any differences. 

Apart from the flux density scale, the flux density 
of the MAGMO data also shows an apparent turnoyer 
around 1.7 GHz. An Incorrect spectral Index In the 
model of PKS 0823-500 could produce such an effect. 
On the other hand, the turnoyer could be real, but It 
would haye to be a yery short-llyed phenomenon, since 
we did not see It In our obseryations in December 2012, 
not eyen a year after the MAGMO ones. The uncon¬ 
firmed result makes It difficult to ascertain the reliability 
of the possible turnoyer. 

In any case, the global flattening of the flux density as 
a function of frequency, which we see in our obseryations 
in a timescale of 1-2 years (Fig. [T]) is a robust result, 
since the same slope is seen in three different frequency 
ranges, obseryed independently. 
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Figure 1. Spectral energy distribution and flux density time-variation of 115103 at radio wavelengths. Our data are represented by solid 
lines: in black the October 2010 data (including filled circles for data at 43-45 GHz); green, August 2011; blue, October 2012 and orange, 
December 2012. The symbols represent archival and literature data as listed in Table 2. Flux density uncertainties in our data are ~ 10%. 
The three points at 4.8 GHz from archives and literature have been slightly displaced in frequency in this figure for clarity, but they were 
taken at the same frequency. 


4. DISCUSSION 

Several scenarios are possible for the origin and vari¬ 
ations of the thermal and non-thermal radio continuum 
emission detected at 115103. We describe in the follow¬ 
ing the different hypothesis and their implication in the 
global context of stellar evolution. 

4.1. Thermal emission from ionized gas 

The weak emission detected on the radio continuum 
map with spectral indices a > —0.1 (Fig. [2|) is com¬ 
patible with the typical free-free emission from ionized 
gas in PNe. It could represent photoionized gas extend- 
ing toward the north east lobe of the mid-IR nebulosity 
(iLagadec et al.ll2qilD, in the direction of the water maser 
let (|G6mez et al.l 120151) . The overall emission in the 
source, however, is strongly dominated by non-thermal 
emission. 

4.2. The origin of non-thermal emission and its 
variation 

A very striking characteristic of the radio emission in 
115103 is that it mostly shows negative spectral indices 
(see Fig. [2]), which indicates that it is dominated by non- 
thermal processes. This emission should arise from elec¬ 
trons moving in a magnetized medium. In particular, the 
spectral index a = —0.54 ± 0.08 found between 8 and 25 
GHz in 2010-2011 is consistent with either synchrotron 
emission from relativistic electrons, or gyrosynchrotron 
from mildly relativistic ones. The spectral index both in 


the case of synchrotron or gyrosynchrotron emission de¬ 
pends on the power-law energy distribution of electrons 
{N{E) oc E~^). Therefore, a change in spectral index, as 
observed between 2010 and 2012, could indicate a time 
variation in this energy distribution that can be used 
to discern which processes are at work. Alternatively, it 
could be the result of absorption or suppression processes 
that are stronger on radiation at lower frequencies. We 
will discuss these possibilities in the following. 

4.2.1. Synchrotron emission and Coulomb collisions 

The spectral index a = —0.54 ± 0.08 is similar to 
that found in synchrotron emissio n in galaxies with ac¬ 
tive nuclei ([Bridle &: PerlevI [19841) . In the case of pure 
synchrotron emission, the spectral index can be related 
to the energy distribution of electrons, in the form a = 
(1 — /3)/2. In our case, the observed a would imply a 
value of /3 = 2.08. 

In the energy conditions of PNe, mass-loss pro¬ 
cesses can produce collisionless, non-relativistic shocks 
and electrons may undergo d iffusive shock acceleration 
(jDrurvlll983l: IAchterber^l200r)() up to very high velocities 
(this process is also called “first-order Fermi accelera¬ 
tion”). In this case, the exponent of the energy distribu¬ 
tion of electrons is given by 


13 = 


r -F 2 
r — 1 


( 1 ) 


where r is the compression ratio of the gas at the shock. 
The maximum value of r in an ideal, non-relativistic gas 
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Figure 2. Top: Map of radio continuum emission integrated be¬ 
tween 8 and 10 GHz. Contour levels are 2'^ (with n from 0 to 
8) times the 3a noise level {a = 1.8 x 10“^ Jy bearn”^). The 
beam size is 1.4'' x 1.1", (position angle = —3°). The thick con¬ 
tour (1.38 X 10“^ Jy beam“^) nearly corresponds to the half-power 
level. The maximum of the image is 2.73 x 10“^ Jy bearn”^. Bot¬ 
tom: Spectral index map obtained with the radio continuum emis¬ 
sion at 8-10 GHz and at 4.5-6.5 GHz, both convolved to a common 
angular resolution of 2.25". Black contours correspond to the spec¬ 
tral indices of —0.1, delimiting the areas that may correspond to 
thermal emission. The white contour is the half-power level of the 
radio continuum map in the top panel. 


is ~ 4. The value /? = 2.08 we obtained indicates that r 
is close to this maximum value. 

In the case of the data taken on October 2012, the 
spectral index of a = —0.28 ± 0.08 found between 4.5 
and 25 GHz would imply an index for the energy distri¬ 
bution of electrons, /3 ~ 1.56 if the emission were due 
to synchrotron emission. However, since the maximum 
compression ratio in collisionless shocks is r = 4, diffusive 
shock acceleration cannot explain indices /3 < 2. There¬ 
fore, the spectral index in the pure synchrotron regime 
should always be a < —0.5, while the acceleration is tak¬ 
ing place. 

However, the synchrotron emission could be the result 
of a shock that took place during a short time and elec¬ 
trons are not being accelerated any longer. In that case, 
the observed spectral flattening reflects the evolution of 
electron energy after the shock e nded. We note that 
115103 is a PN (|G 6 mez et akllMsIl . and it has developed 
a photoionized region around it. When the relativistic 


electrons move in a dense plasma, they lose energy due 
to Goulomb collisions against ions. The lifetime of elec¬ 
trons of_Jov^energy_Js_shurterunder Coulomb losses 
(e.g., iRephaelilllOTol : [^razinlll999ll : 


icou(7) ^ 7 1-2 X 10 




1.0 


ln(7/ne 


75 


s, 


( 2 ) 

were Ue is the electron density of the plasma expressed in 
cm“^, and 7 is the Lorentz factor of the relativistic elec¬ 
trons. The result is a progressive hardening of the energy 
distribution, with less electrons of low energy. This will 
lead to a flattening of the synchrotron spectrum. Spec¬ 
tral flattening due to Coulomb collis ions is also observed 
in the radio spectrum of solar flares (IMelnikov fc Maeunl 
[l^lLee fc C^IMlOll . 


The lifetimes of electrons under Coulomb collisions are 
consistent with the observed flattening of 115013 over 
time. In particular, assuming a mag netic field of ~ 5 
mC, as seen in so me po s t-AG B stars (jBains et al.ll2003L 
1200411 . according to IShul (|1991ll . the Larmor frequency of 
electrons 


eB 

27r rn-eC 


(3) 


(where e and me are the electron charge and mass, re¬ 
spectively, and B is the magnetic field strength) would be 
~ 14 kHz. Thus, the synchrotron emission at cs 1 GHz 
would typically be due to electrons with 7 = ^JvJvl = 
267. 

Equation[2] indicates that the lifetime of these electrons 
under Coulomb collis ions and assuming electron density 
of 2 X 10^ cm“^ (e.g. iTafova et al.ll20Cil^ would be ~ 40 
yr. This is a reasonable value considering the timescales 
of variations at low frequency shown in Fig. [TJ Under 
this interpretation, the observed non-thermal emission in 
115103 would be due to synchrotron emission from elec¬ 
trons accelerated during an explosive event that occurred 
a few years ago, and the emission would now be decaying 
due to Coulomb losses. 

We note that a recent explosive event would also ex¬ 
plain t he velocity gradient of the maser emission in this 
source (jComez et al.11201^ . which reinforces this inter¬ 
pretation. 

We also note that the timescale of relativistic electrons 
will shorten with increasing plasma density (equation [5]). 
The change in spectral index between 2010 and 2012 
seems particularly sudden compared with the evolution 
in previous years. This is suggestive of an increasing den¬ 
sity of ionized gas, consistent with 115103 being a very 
young PN that is developing its photoionized region. 


4.2.2. Gyrosynchrotron emission from mildly relativistic 

electrons 

An alternative explanation to the observed change in 
the SED of 115103 would be that we witness the variation 
in gyrosynchrotron emission from mildly relativistic elec- 
trons. This radiative process is seen in active stars (e.g., 
lDulklll98'5l : lCunn et al.lll993 : lPestalozzi et all2000[ l. The 
flux density of these stars roughly follows a power-law of 
the form 


5, OclO-°-52/5(sin0)-°-43+O-65/3 



(4) 
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where the relation between the spectral index and the 
distribution of electrons is given by a ~ 1.22 — 0.90/3 
9 is the pitch angle between the electron 
trajectory and the magnetic field, and vi, ~ 2.8 x 10®B 
is the electron-cyclotron frequency. 

Equation (4) does not take into account any effect 
from the thermal plasma, and gyrosynchrotron emis¬ 
sion will in general deviate from a s i mple power law 
(see, e.g. Fleishman fc Kuznetsovlll2010l: iKuznetsov et al.l 
1201111 . However, it serves to illustrate the feasibility of 
this mechanism in the case of 115103. 

As in the previous hypothesis, the change in spectral 
index would imply a change in the energy distribution of 
electrons, and all spectral indices observed in 115103 are 
in principle possible. Considering that the flux density 
has not changed much at 23 GHz between 2010 and 2013, 
equation |4] implies that magnetic fields of several kG are 
necessary to flatten the radio spectrum while maintaining 
a similar flux density at 23 GHz. These m agnetic field 
stren gths have not been observed in PNe (ILeone et al.l 
l2014fl . which makes this hypothesis highly unreliable. 
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Figure 3. Illustration of the effect of the opacity due to ionized 
gas. Dashed line: linear fit to the data on 2010-2011. Dotted 
line: linear fit to the data on October 2012. Dashed-dotted line: 
Expected emission by applying the opacity of ionized gas (r oc 
to the fit to the 2010-2011 data. Color codes and symbols 
for the data are the same as in FigU 
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4.2.3. Opacity of free thermal electrons and Razm effect 

We have also explored the possibility that the change in 
the spectral index were not due to a variation in the elec¬ 
tron distribution, but to selective absorption/suppression 
of radiation at low frequencies. If the emission on 2010- 
2011 is of gyrosynchrotron/synchrotron nature, this can 
be later suppressed by an increasing amount of ionized 
gas via two possible processes: a) opacity from free 
electrons or bj the Raz i n for Razin-Tsytovich) effect 
(iRazinI 119601 : iTsvtovichl 119511 : iGinzburg fc Svrovatskiil 
119651) . Both processes are more efficient at lower fre¬ 
quencies. 

In the former case (a), the opacity due to free electrons 
has a dependency oc We tested this possibility 

against our data, assuming a screen of ionized material 
between the observer and the source of non-thermal emis¬ 
sion, and applying an optical depth Ti, = to a SED 

with spectral index —0.54, varying the constant a until 
we reached a final spectral index of ~ —0.28 between 8 
and 23 GHz. However, the resulting screen of ionized 
gas would have such a high opacity at lower frequencies 
that we should not detect any emission at ~ 1 — 2 GHz 
(Fig. [3]). Moreover, an optically thick screen of free elec¬ 
trons would produce thermal free-free emission with a 
flux density proportional to where H is the solid 

angle subtended by the screen. Given the frequency de¬ 
pendence of free-free emission, it cannot compensate the 
deficit of emission at 1-2 GHz since, in that case, ther¬ 
mal emission would dominate at even higher frequencies, 
with a positive spectral index that is not observed. 

In the case of Razin effect (b), the shape of the spec¬ 
tral energy distribution of synchrotron emission within a 
plasma takes the form 

P(e, /3) = jp ^(/3-3)/2^^ ^5) 

(function P in eauation[51in lPacholczvklll970L but with 
some changes in the nomenclature of variables), where 

/=(! + -) (6) 


The function F is defined as 

nOC 

F{y) = y K^/^{z)dz (8) 

Jy 

where ATs/s is the Bessel function of order 5/3. 

The key parameter is the ratio n^jiB sin 0), which con¬ 
trols the change in spectral slope. Assuming sind ~ 1, 
the critical frequency below which the synchrotron emis¬ 
sion is mostly suppressed is ~ 2^n^/B. Applying Equa¬ 
tion [5] to a SED with spectral index —0.54, to reach a 
final spectral index of ~ —0.28 between 8 and 23 GHz, 
would also imply the complete suppression of emission 
at low frequency, as in the case of plasma opacity (a), 
ruling out these two hypothesis. 

4.2.4. Combination of thermal and non-thermal emission 

A possible explanation for the observ ed change in spec- 
tral index is the behavior predicted bv iFleishmanI (120061 1 
for non-thermal emission within dense plasmas when the 
magnetic field has both an ordered and a random com¬ 
ponent. 

We note that the spectral index a ~ —0.28 ob- 
served in 2012 is simi lar to that in the Grab nebula 
( Bietenholz_et.jl[J]^97|) and other pulsar wind nebulae 
(|Gaensler fc Sland 120061 1. In these cases, the emission 
can be int erpreted as “jitter” (or “diffus i ve synchrotron” 1 
radiation iTootvein fc Fleishmanl[l987t lMedvede'^l2000t 
iFleishman fc Bietenhold 1200711 . i.e.. electrons moving in 
a random magnetic field, as opposed to the ordered field 
for usual synchrotron emission. 

If we consider 115103 as having both an ordered and a 
random component, where the ordered component dom¬ 
inates, relativistic electrons moving along the magnetic 
field lines would mainly show a normal synchrotron emis¬ 
sion, with a relatively small contribution of jitter emis¬ 
sion. However, if this magnetic field is embedded in a 
plasma, the Razin effect would suppress the synchrotron 
component. The suppression of jitter emission depends 
much more slow ly on plasma density (equation 40 in 
iFleishm^ 1200611 . and therefore, if this density is suffi¬ 
ciently high, the jitter emission will eventually dominate 
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the entire spectrum even though the magnetic field and 
the electron energetics remained unchanged. The Razin 
effect is actually stronger in the presence of dense plas¬ 
mas and relatively weak magnetic fields, so we expect it 
to be more readily detectable under the physical condi¬ 
tions of PNe, over those in other synchrotron-emitting 
objects, such as radio galaxies or supernovae. 

In this scenario, and for the case of 115103, until 2010, 
while the photoionization of the PN was still low, it would 
show a mixture of synchrotron and jitter emission (con¬ 
sistent with a magnetic field with an ordered and a tur¬ 
bulent component). The increasing plasma density of 
the newly formed ionized region of the PN would then 
suppress the synchrotron component due to Razin effect. 
The final result would be pure jitter emission in the 2012 
data. If the change in the spectral slope is due to a 
sudden change in the dominant radio continuum emis¬ 
sion mechanism, it would suggest that we are witnessing 
the onset of the photoionization in a PN. However, given 
that we do not have a direct determination of the mag¬ 
netic field strength, the plasma density, and the size of 
the non-thermal emitting region, it is not possible at this 
point to determine whether this interpretation is correct. 

It is interesting to compare the flattening in the spec¬ 
tral energy distribution of the radio continuum emis¬ 
sion, with th e flattening seen in the post -AGB star IRAS 
15445—5449 (|Perez-Sanchez et al.ll201§) . whose nature is 
clearly different from the one we see in 115103. In IRAS 
15445—5449, the change in the slope between 2.5 and 22 
GHz is due to an increase of flux density at high frequen¬ 
cies. This was interpreted as an increasingly higher con¬ 
tribution of free-free emission which could arise from a 
shock-ionized wind. It is known that a highly collimated 
ionized wind ca n produce significant free-free emission 
(jRevnoldiJ 119861 1. and a stronger emission at high fre¬ 
quencies could be the consequence of a sudden increase 
in its mass-loss rate. On the contrary, the flattening we 
observe in 115103 is due to a decrease of flux density at 
low frequencies. The flux density distribution seems to 
go down, pivoting around u ~ 23 GHz, where it does not 
change much. This cannot be explained by an increase of 
free-free emission, but strongly indicates a suppression of 
emission at low freque ncies. Other post- AGB stars show 
non-thermal emission (jBains et al.ll2009ll , but there is no 
information on their time evolution, so a comparison is 
not possible at this stage. 

Flat radio spectra with ind ices a > —0.5 have also been 
found in Wolf-Rayet stars (jOhaoman et al.l [19990 . and 
explained as a combination of thermal and non-thermal 
emission in colliding wind binaries. Under these models, 
the non-thermal component has a < —0.5, and the flat¬ 
tening would be attributed to the thermal mechanism. 
However, as mentioned above, we cannot explain the 
progressive flattening of the spectrum in 115103 with an 
increase of thermal emission. 

4.3. Processes taking place during the birth of a 
planetary nebula 

115103 is the first known PN in which synchrotron 
emission could be detected, since the free-free emission 
is so weak that it has not masked the non-thermal emis¬ 
sion yet. This non-thermal emission due to shocks could 
share a common origin with that seen in some post -AGB 
stars (|Bains et al.ll2009HPerez-Sanchez et al.ll2013ll . 


However, in the case of 115103, if we assume the 
Goulomb losses as the explanation for the variability in 
the spectral index (Sec. 14.2.IL we would be observing 
the results of a jet launched close to the moment when 
the post-AGB becomes a PN. Theoretical predictions af¬ 
firm that jets can be launched when magnetic field lines 
are twisted due to the different rotation speeds of the 
envelope and the stellar core of an evolved star, and the 
energy is released in the form of a collimated ejection 
(jMatt et al.ll200fill . The conditions for this magnetic ex¬ 
plosion will be favoured when a star enters the PN phase, 
as is the case of 115103, since the magnetic field lines 
away from the star will tend to be locked by the ion¬ 
ized gas to the slower rotational motions of the outer 
envelope. Moreover, recent studies show that magnetic 
fields are stronger in older post-AGB stars and might 
be related to t he origin of non-spher ical morphologies in 
these objects (jGonidakis et al.ir2014fl . 

The changes in the spectral index could also be pro¬ 
duced by the Razin effect (Sec. I4.2.4L where the be¬ 
ginning of photoionization could start new physical pro¬ 
cesses that affect the behavior of the radio continuum 
emission. 

Previous studies of radio c ontinuum in very young PNe 
showed that in severa l cases (|Cerrigone et al.l[2008L l20lil : 
iGasassus et al.l[200^ the observed spectral index could 
not be explained by the classical optically thin free-free 
emission and the deviations from the expected a were 
in some cases attributed to calibration uncertainties or 
other instrumental effects. Radio continuum variability 
as well as negative spectral indic es were also found in 
young PNe (iCerrigone et al.l[2011[l . but limited time and 
frequency coverage prevented a detailed identification of 
the processes taking place. In fact, it is possible that 
some of these young PN may be undergoing the same 
processes as we are witnessing in 115103, but it has only 
been possible to identify them in this study, thanks to the 
implementation of new broad-band receivers and corre¬ 
lators in radio telescopes such as ATGA or the Jansky 
Very Large Array. 

4.4. Predictions on source characteristics and their 
future evolution 

A monitoring of the radio spectrum of 115103 would 
shed light on the processes taking place in this unique 
object. As shown above, possible origins of the variation 
seen in the spectral index are changes in synchrotron 
emission due to Goulomb collisions (Sec. 14.2.111 . or se¬ 
lective suppression of part of the emission due to Razin 
effect (Sec. 14.2.41) . However, we expect that these pro¬ 
cesses would lead to different future evolutions in the 
short term. In the case of changes in the energy distri¬ 
bution of electrons due to Coulomb collisions, we expect 
that the spectral flattening will continue. Eventually, a 
low-frequency cutoff will arise in the spectrum, reflect¬ 
ing a complete suppression of electrons below an energy 
threshold. On the other hand, if the physical process is 
Razin suppression of synchrotron emission, leaving only 
jitter emission, we do not expect a significant further 
flattening of the emission, since the dominant radiative 
process after 2012 should not be immediately affected by 
the Razin effect. 

As seen in Fig. [2l the possible contribution of free- 
free emission from ionized material is still much weaker 
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than the non-thermal emission in this source. In the long 
term, we predict that free-free radiation will significantly 
increase in the future, as photoionization proceeds. Its 
more immediate effects would be an increase of flux den¬ 
sities at high frequencies due to free-free emission, and 
a decrease at low frequencies due to the absorption of 
non-thermal emission by the plasma opacity. At later 
times, free-free emission will eventually dominate the ra¬ 
dio spectrum of the source, giving rise to a ther mal SED 
similar to other young PNe (|G6mez et al.l[2f)?)^ . with no 
trace of the non-thermal emission that we observe now. 

5. CONCLUSIONS 

We have observed radio continuum emission between 1 
and 46 GHz towards the young PN IRAS 15103—5754, in 
different epochs between 2010 and 2012, using the Aus¬ 
tralia Telescope Compact Array. Archive radio contin¬ 
uum data from 1991 to 2012 have been used to complete 
the information on this source. Our main conclusions are 
as follow: 

• Radio continuum emission shows a non-thermal 
spectrum, which is not typical of PNe. 

• The flux density seems to have decreased with time, 
especially at the lowest frequencies. In particular 
we see a flattening of the spectrum, with a spec¬ 
tral index a ~ —0.54 in 2010-2011, which changed 
to a ~ —0.28 in 2012. While the radio spectrum 
observed in 2010-2011 can be explained as arising 
from synchrotron emission with electrons acceler¬ 
ated in a shock, the one in 2012 is not compatible 
with this process alone. 

• The variation of the radio spectrum in 115103 can 
be attributed to two effects of the onset of pho¬ 
toionization in this source: 

— The loss of energy of the synchrotron-emitting 
electrons due to Coulomb collisions within the 
ionized region of the PN. The estimated life¬ 
time of the electrons is 18 yr, consistent with 
the timescale of variation of the flux density 
in the archival data. The onset of the pho¬ 
toionization could produce the sudden change 
visible in the spectral index between 2010 and 
2012 . 

— The electrons moving in a magnetic field with 
both an ordered and a random component, 
which would give rise to a mixture of syn¬ 
chrotron and ’’jitter” emission. The Razin ef¬ 
fect would suppress the synchrotron compo¬ 
nent as photoionization progresses. Only the 
jitter component with a flatter spectral index, 
will remain. 

• We postulate that the flattening of non-thermal 
spectra could be a signature of the onset of the 
PN phase. 
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